Spin- Wave Description of Haldane-gap antiferromagnets 
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Modifying the conventional antiferromagnetic spin-wave theory which is plagued by the difficulty 
of the zero-field sublattice magnetizations diverging in one dimension, we describe magnetic prop- 
erties of Haldane-gap antiferromagnets. The modified spin waves, constituting a grand canonical 
bosonic ensemble so as to recover the sublattice symmetry, not only depict well the ground-state 
correlations but also give useful information on the finite-temperature properties. 



PACS numbers; 75.10.Jm, 05.30.Jp, 75.40.Mg 

Haldane's conjecture [1,2] that one-dimensional 
Heisenberg antiferromagnets should exhibit qualitatively 
different low-energy structures based on whether the con- 
stituent spins are integral or fractional sparked renewed 
interest in the field of quantum magnetism. An energy 
gap immediately above the ground state was indeed ob- 
served in a quasi-one-dimensional Heisenberg antiferro- 
magnet Ni(C2H8N2)2N02(C104) [3] and a rigorous ex- 
ample of such a massive phase was also found out [4,5]. 
The energy gaps in magnetic excitation spectra, that is, 
spin gaps, are now one of the most attractive and im- 
portant topics. In the context of theoretical progress, we 
may be reminded of quantized plateaux in the ground- 
state magnetization curves [6] , a dramatic crossover from 
one- to two-dimensional quantum antiferromagnets [7], 
and an antiferromagnetic excitation gap accompanied by 
ferromagnetic background [8-11]. From the experimen- 
tal point of view, metal oxides such as spin-Peierls com- 
pounds Cui_^M^Ge03 (M = Zn,Mg) [12,13], Haldane- 
gap antiferromagnets i?2BaNi05 (i? = rare earth) [14] 
and ladder materials Sr„_iCu„+i02Ti {n — 3, 5,7, •••) 
[15] have significantly contributed to systematic investi- 
gations of the mechanism of gap formation. 

In studying Haldane-gap antiferromagnets, numeri- 
cal tools such as quantum Monte Carlo and density- 
matrix renormalization-group methods are indeed use- 
ful, but analytic approaches still play an important role. 
The nonlinear tr-model technique [1,2,16-19] is powerful 
enough to investigate low-energy structures with partic- 
ular emphasis on the competition between massive and 
massless phases. It is the valence-bond-solid description 
[4,5] of integer-spin chains that enables us to readily un- 
derstand the hidden order [20] inherent in the Haldane 
massive phase. However, these methods are not so use- 
ful for exploring thermal properties as they are in in- 
vestigating the ground-state properties. Although we 
may consider applying the spin- wave scheme to Haldane- 
gap antiferromagnets in such circumstances, the diffi- 
culty of the zero-field sublattice magnetizations diverg- 
ing in one dimension has been a problem in the conven- 
tional spin-wave theory [21,22] for years. Therefore, it 
was a major breakthrough that spin waves, being con- 
strained to keep zero magnetization, succeeded in pre- 
cisely describing the low-temperature thermodynamics of 



the one-dimensional Heisenberg ferromagnet [23]. This 
new spin- wave scheme, which is now referred to as the 
modified spin-wave theory, was further applied to quan- 
tum antiferromagnets and ferrimagnets. The modified 
spin-wave scheme is highly successful for extensive ferri- 
magnets [24-29] and still applies well for two-dimensional 
antiferromagnets [30-34] . As for one-dimensional antifer- 
romagnets, there exists a pioneering argument [35], but 
it looks unsatisfactory for interpreting experimental and 
numerical observations. 

In this article, we aim to do away with our vague but 
persistent impression that the spin-wave scheme hardly 
works for one-dimensional quantum antiferromagnets. 
Spin-2 Haldane-gap antiferromagnets [36,37], as well as 
those of spin 1, have recently been synthesized and more 
explorations into novel quantum phenomena in one di- 
mension are expected in the future. In such circum- 
stances, we make our first attempt to construct a modified 
spin-wave theory for Haldane-gap antiferromagnets. 

We consider integer-spin antiferromagnetic Heisenberg 
chains 

L 

■H = jJ^Sj ■ Sj+i; S^j=S{S + l). (1) 

We define bosonic operators for the spin deviation in each 
sublattice via 
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The Fourier-transformed operators are introduced as 
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6fc = ^^e-^-("+i/4)6„, 

n=l 



where twice the lattice constant 2a is set equal to unity 
and therefore k = 2Tm/N (n = 0, 1, • • • , N-l; N = L/2). 
Through the Bogoliubov transformation 



a\ = a[.cosh0fe — /3fcsinh6'/j . 
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we obtain 

n = E^ + Ex + Eo + ni + no + o{s-^) , 

where 

= -2NJS'^ , 
El = -ANJS{r - A) , 
Eo = -2NJ{r - Af , 

Hi = J^[wi(A;) (alak + PlPk) 
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7i(A;) — 2S ( cos — cosh2^?/j. — sinh20fe 



70 (fc) = 2(r - ^) ( cos - cosh26'fe - sinh26»fe ) . 
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The naivest diagonahzation of the Hamiltonian (5), 
whether up to 0{S^) or up to 0(5'^), results in di- 
verging sublattice magnetizations even at zero temper- 
ature. Although cq. (2) assumes that a„ < 25* and 
bnbn < 25, the conventional spin-wave theory cannot 
reasonably control the boson numbers. Then we consider 
introducing a grand canonical constraint to the noncom- 
pact Hamiltonian (5). 

Isotropic magnets should lie in the state of zero mag- 
netization Sj = and the minimization of the free 
energy under such a condition indeed yields an excellent 
description of the low-temperature thermodynamics of 
ferromagnets [23]. However, in the cases with antiferro- 
magnetic exchange interactions, the zero-magnetization 
constraint, claiming that X)n('^n'^" ^b\hn) — 0, still fails 
to overcome the divergence of the numbers of the sub- 
lattice bosons. Hence we minimize the free energy con- 
straining the sublattice magnetizations to be zero: 
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= Y,bibn = SN. 



(11) 



Within the conventional spin-wave theory, spins on one 
sublattice point predominantly up, while those on the 
other predominantly down. This condition (11) restores 
the sublattice symmetry. In order to enforce the con- 
straint (11), we first introduce a Lagrange multiplier and 
diagonalize 



n=n + 2JXj2iaW + blbk) . 



(12) 



Then the ground-state energy and the dispersion relation 

are obtained as 

Eg^ E^+Ei; Ei^Ei+ ANJXA , (13) 
uj{k) = uji (k) ; LOi (k) = loi (k) + 2Acosh26'fe , (14) 

within the linear modified spin-wave scheme and as 



Eg = i^N + El + Eq , 
w(fc) — uji{k) + LOo{k) 



(15) 
(16) 



within the up-to-0(5'°) interacting modified spin-wave 
scheme. For eqs. (13) and (14), 9k is given by 71 (A;) — 
2Asinh2^?fc = 71 (fc) — 0, whereas for eqs. (15) and 
(16), 9k may be determined in two ways. One idea is 
the perturbational treatment of Ho, which is referred 
to as the perturbational interacting modified spin-wave 
scheme, where 9k is still given by 71 (fc) = and the 
0{S^) quantum correction is the 0(S'^)-eigenstate av- 
erage of Hq. The other is the full diagonahzation of 
Til + Ti.0, which is referred to as the fuU-diagonalization 
interacting modified spin-wave scheme, where 9k is given 
by 7i(fc) + 7o(fc) = 0. Once 6k is given, we calculate the 
free energy and obtain the optimum thermal distribution 
functions as 



(alak) = iPipk) 
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(17) 



where A is self-consistently determined by the condition 
^(2nfe-|-l)cosh26»fc = (25-1-1)^^. (18) 

k 

First, let us evaluate the ground-state energy. We com- 
pare the modified spin-wave calculations with the highly 
accurate quantum Monte Carlo estimates [38] in Table 
I. The modified spin-wave findings are generally in good 
agreement with the quantum Monte Carlo results. The 
interacting modified spin waves describe the ground-state 
correlations much better than the linear ones and their 
description becomes increasingly refined with increasing 
S. The interacting modified spin-wave findings miss the 
correct value by 0.5 percent for 5 = 1 and by only 0.008 
percent for 5 = 3. 

Secondly, we consider the Haldane gap w(afc = tt) = 
Z\(T). Table II shows that the modified spin- wave 
scheme, in contrast with the conventional spin-wave the- 
ory, succeeds in generating the gap but considerably un- 
derestimates it. Unavailability of the absolute energy 
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level structure is an inevitable consequence of our em- 
ploying the effective Hamiltonian (12). Then we ex- 
amine the present scheme by scaling A(T) to its zero- 
temperature value A(T = 0) = Aq. Such an argu- 
ment is quite usual with field-theoretical investigations, 
which do not lead to an estimate of the normalization 
factor but derive finite-temperature expressions involv- 
ing only ratios such as k^T/A^. In Fig. 1, we com- 
pare the modified spin-wave calculations with numerical 
[39], field-theoretical [40] and experimental [41] findings. 
Now the superiority of the modified spin- wave scheme is 
clear at a glance. Of all the theoretical tools, the fuU- 
diagonalization interacting modified spin-wave approach 
is the most successful to reproduce the observed upward 
behavior of A{T) with increasing temperature [41-43]. 
The nonlinear cr-model treatment [40,44] is justified well 
in the low-temperature region k^T <C Af), while the 
maximum-entropy technique [39] works less with increas- 
ing temperature. Nickel compounds such as Y2BaNi05 
[41] and Ni(C2H8N2)2N02(C104) [42] are good candi- 
dates for spin-1 Haldane-gap antiferromagnets, but mag- 
netic anisotropy and interchain interaction are not neg- 
ligible there. They split the lowest excitation gap of the 
isotropic chain (1) into several levels and smear the in- 
trinsic behavior of the ideal integer-spin chains. Consid- 
ering such practical factors, the present theory satisfac- 
torily interprets the observations. 

Lastly, we show the modified spin-wave calculations of 
the magnetic susceptibility in Fig. 2. Considering the 
total breakdown of the conventional spin-wave theory in 
one-dimensional thermodynamic calculations, the modi- 
fied spin-wave achievement is highly successful. All the 
calculations are guaranteed to reproduce the paramag- 
netic susceptibility x/iff^Me = + l)/3fcBr at high 
temperatures. Since the interacting modified spin-wave 
scheme gives better estimates of the gap than the linear 
modified spin- wave one (Table I), it is somewhat bet- 
ter at describing the low-temperature behavior. With 
increasing S*, the activation-type initial behavior is sup- 
pressed and the antiferromagnetic peak is broadened. 
The S = 1 modified spin-wave calculations are in fine 
agreement with the quantum Monte Carlo findings over 
a wide temperature range, while those for S = 2 look 
somewhat poorer at intermediate temperatures. It may 
be closely related to the fact that the modified spin- wave 
estimates of the 5 = 2 gap are worse than those of the 
S = 1 gap (Table II). However, Table II suggests that 
the validity of the modified spin-wave scheme for excita- 
tions significantly improves with increasing S*, possibly in 
a staggered way. Thermodynamic calculations, whether 
by quantum Monte Carlo or density-matrix renormaliza- 
tion group, for systems with large degrees of freedom are 
less feasible numerically, in particular, at low tempera- 
tures. The present scheme has the advantage of saving 
time and computational resources. 

We have demonstrated the applicability of the new 
spin-wave scheme to Haldane-gap antiferromagnets. This 
is the first comprehensive attempt to describe one- 



dimensional spin-gapped antiferromagnets in terms of 
spin waves. The modified spin waves are free from their 
thermal as well as quantum divergence and can there- 
fore microscopically interpret various magnetic proper- 
ties. Besides the magnetic susceptibility, the spin cor- 
relation function and the nuclear spin-lattice relaxation 
time can be revealed. Unfortunately, we have less in- 
formation on the bare energy spectrum, because we get 
rid of the quantum divergence at the cost of the original 
Hamiltonian. However, with the zero-temperature spin 
gap, which is readily and precisely available through nu- 
merical calculations within a canonical ensemble [38] , we 
can still argue the energy structure quantitatively. 
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FIG. 1. Temperature dependences of the spin-1 Haldane 
gap calculated by the linear (LMSW), perturbational interact- 
ing (PIMSW) and full-diagonalization interacting (FDIMSW) 
modified spin-wave schemes, a quantum Monte Carlo (QMC) 
method combined with the maximum-entropy technique [39] 
and the nonfinear cr-model (NSM) approach [40]. Inelastic 
neutron-scattering measurements (Exp.) on Y2BaNi05 [41] 
are also shown for reference. 
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FIG. 2. The linear (LMSW), perturbational interacting 
(PIMSW) and full-diagonalization interacting (FDIMSW) 
modified spin-wave calculations of the zero-field magnetic sus- 
ceptibility as a function of temperature compared with quan- 
tum Monte Carlo (QMC) estimates. 
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Haldanc-gap antifcrromagnets in a magnetic field pro- 
vide further interesting issues. With the increase of an 
applied field, the gap is reduced and the ground state is 
mixed increasingly with the first excited state. Indeed the 
nuclear spin-lattice relaxation is accelerated in the vicin- 
ity of the critical field, but its overall behavior as a func- 
tion of an applied field and temperature [45] is far from 
understandable at a glance. We may expect the present 
new scheme to open the way for the total understanding 
of low-dimensional spin-gapped antiferromagnets. 
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TABLE I. The Hnear (LMSW), perturbational interact- 
ing (PIMSW) and full-diagonalization interacting (FDIMSW) 
modified spin-wave calculations of the ground-state energy 
per site compared with quantum Monte Carlo (QMC) esti- 
mates [38]. 





5=1 


5 = 2 


5 = 3 


LMSW 


-1.361879 


-4.726749 


-10.0901 


PIMSW 


-1.394853 


-4.759760 


-10.1231 


FDIMSW 


-1.394617 


-4.759759 


-10.1231 


QMC 


-1.401481(4) 


-4.761249(6) 


-10.1239(1) 



TABLE II. The linear (LMSW), perturbational interact- 
ing (PIMSW) and full-diagonalization interacting (FDIMSW) 
modified spin-wave calculations of the lowest excitation gap 
A{T = 0) = Ao compared with quantum Monte Carlo (QMC) 
estimates [38]. 





5=1 


5 = 2 


5 = 3 


LMSW 


0.07200 


0.00626 


0.00279 


PIMSW 


0.07853 


0.00655 


0.00287 


FDIMSW 


0.08507 


0.00683 


0.00295 


QMC 


0.41048(6) 


0.08917(4) 


0.01002(3) 
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